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ABSTRACT 

Aims. There are large differences in the proposed progenitor models for the Cas A SNR. One of these differences is the presence or 
absence of a Wolf-Rayet (WR) phase of the progenitor star. The mass loss history of the progenitor star strongly affects the shape of 
the Supernova remnant (SNR). In this paper we investigate whether the progenitor star of Cas A had a WR phase or not and how long 
it may have lasted. 

Methods. We performed two-dimensional multi-species hydrodynamical simulations of the CSM around the progenitor star for several 
WR life times, each followed by the interaction of the supernova ejecta with the CSM. We then looked at the influence of the length 
of the WR phase and compared the results of the simulations with the observations of Cas A. 

Results. The difference in the structure of the CSM, for models with different WR life times, has a strong impact on the resulting 
SNR. With an increasing WR life time the reverse shock velocity of the SNR decreases and the range of observed velocities in the 
shocked material increases. Furthermore, if a WR phase occurs, the remainders of the WR shell will be visible in the resulting SNR. 
Conclusions. Comparing our results with the observations suggests that the progenitor star of Cas A did not have a WR phase. We 
also find that the quasi-stationary flocculi (QSF) in Cas A are not consistent with the clumps from a WR shell that have been shocked 
and accelerated by the interaction with the SN ejecta. We can also conclude that for a SN explosion taking place in a CSM that is 
shaped by the wind during a short (< 15000 yr) WR phase, the clumps from the WR shell will be visible inside the SNR. 

Key words. Hydrodynamics - stars: winds, outflows - stars: supemovae: general - stars: supemovae: individual: Cas A 



■ 1 . Introduction limit of the mass of a star at the start of the m ain sequence (MS) 
^ „ , , „ „ ,oxtt^n ■ • in the Galaxy that can b ecome WR stars (Meynet & Maeder 
One of the challenges of supernova remnant (SNR) research is 2QQ ^ | Ma ssev et al.ll200(ih . However, it is not clear whether the 
^ , to link the properties of the remnant with the properties and the bw ejecta mass ig due tQ h mass 1qss of a ^ massiye 
evolution of the progenitor star. Ideally one would hope to learn star ^ was Qn ife or just entered the WR phase, or whether 
about the late stages of stellar evolution from the properties of mass lo ss was induced by a common envelope phase in a binary 
the SNR In that respect one of the best-studied SNRs is the system dYovmg et aLl lM. Als o the presence of a jet-counter- 
bright Galactic remnant Cas A. jet syfjtem jn Cas A ^^2004 iHwang et al]l200l puts strong 
For a long time Cas A has been thought to be the rem- constraints on the duration of a possible WR phase Ischure et all 
nant of a very massive star, possibly that of an exploding Wolf- rj 008) . Furthermore, the high density of the shocked circumstel- 
Rayet (WR) star (e.g.| Fesen et alJ | 1987 fc| Woosley et al|1993|) . In lar material (CSM) indicates that the shock wave is cuirently 
that case one would expect the supernova to be of Type Ib/c. moying throug h the red super g ia nt (RSG) wind of the progeni- 
However, recently IKrause et alj d2008l) obtained a spectrum of tor dchevalier & Oishi 2003|) 
the supernova by observing its light echo, identifying it as a Type 

lib explosion. This indicates that the progenitor lost most, but Although recent studies dSchure et alj l2Q08h IKrause et all 

not all of its hydrogen envelope. I2008I) seem to suggest an extended progenitor star, the pres- 

X-ray observations indicate that the ejec ta mass was only ence of a WR shell, i.e. the shell of red supergiant (RSG) wind 

2 - 4 M n with about 1 - 2 M G of oxygen dVink et all 119961; material swept up by the fast WR wind, has been invoked to 

IWillingale et al.ll2003l) . The latter suggests an initial mass for the explain the presence of numer ous slow moving (< 500 km/s 

progenitor of approximately 20 M Q , which is close to the lower Ivan den Bergh & Kampe rH985| ) nitrogen rich knots often c alled 

quasi-stationary flocculi (OSF). IGarcfa-Segura et al.l ( 1996) sug- 

Send offprint requests to: B. van Veelen gested that these knots are the remnants of the broken up WR 

* See the online version for color figures. shell. Note that also a binary common envelope phase is likely to 
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be followed by a WR-like mass loss phase ( Young et al. 2006). 
Both the fact that the QSF lie within the boundary of the for- 
ward shock, and the inferred density behind the forward shock 
of Cas A, limits the possi ble duration of a WR phase to ~ 10 4 yr 
dGarcfa-Segura et alJ[l996b . 

Here we report on our investigation of the imprint of the pro- 
genitor's stellar wind evolution, i.e. the existence and duration of 
a WR phase, on the morphology and kinematics of Cas A. For 
the first time we study the hydrodynamical evolution of a SNR in 
2-D using multi-species advection, which allows us to separate 
the location and kinematics of wind material from the supernova 
ejecta. 

2. Method and Assumptions 

2. 1 . The adopted stellar evolution model 

As discussed in Section [T] the progenitor star of Cas A has 
evolved through a RSG phase and perhaps a WR phase before 
exploding. This evolutionary history gives us constraints on the 
mass of the progenitor. Unfortunately, the allowed mass range of 
stars th at have a WR phase during their evolution is not very well 
known. Meynet & Maeder (2005) show that stars with an initial 
mass as low as 20M Q may end up as a WR star, and that those 
stars will have a short WR life time. Alternatively, the progeni- 
tor of Cas A may have gone through a binary common envelope 
phase. However, we treat here the evolution of the progenitor as 
that of a single massive star. The reason for that is partly prac- 
tical; there are no good analytical models of the hydrodynamics 
of a common envelope phase, and a full numerical simulation of 
this is beyond the scope of this study. Moreover, one may expect 
that a common envelope phase results in a non-spherically sym- 
metric outflow. Although there is ev idence for non-spheric ity of 
the supernova ejecta in Cas A (e.g. [Willingale et al. 2002J), the 
outer shock wave is surprisingly circular (IVinkl i2004). suggest- 
ing a more or less spherically symmetric CSM. The CSM struc- 
ture that we assume, i.e. one that is shaped by a spherically sym- 
metric wind of a single massive star, seems valid for Cas A. 

We have chosen an initial mass for the progenitor star of 
20Mq, in accordance with previously suggested ZAMS masses 
dVinkll2004 . The stellar evoluti on model of this 20M Q , non- 
rotating star was calculated by iHirschi et al.l d2004l) . In this 
model the star does not have a WR phase and we thus assumed 
an enhanced RSG mass loss in such a way that during the RSG 
phase the entire hydrogen envelope is lost. It does not matter for 
the evolution of the core if the star loses its envelope entirely or 
only partly in this case, which means that we can assume the 
enhanced mass loss without affecting other parameters of the 
stellar evolution model. An enhanced mass loss rate is not un- 
reasonable given by the uncertainty there currently is in RSG 
mass loss rates (see Ide Jager et al.lll988t Ivan Loonetal.ll200l 
and discussion therein). 

The core mass at the end of the RSG phase is 6M which 
is also the final mass of the progenitor, due to our assumed en- 
hanced mass loss rate. This mass is also roughly consistent with 
the ejecta mass of Cas A of about 2 - 4M Q , taking the pres- 
ence of the neutron star in Cas A into account dTananbauml 1 999t 
Chakraba rtv et al.1 12001). We adopt an ejecta mass of 4M Q , and 
an ejecta kinetic energy of 2 ■ 10 erg dVinkll2004l) . 

The mass loss rate and velocity of the stellar wind are cal- 
culated by combining the observational constraints and the pa- 
rameters of the stellar evolution model. We can derive the total 
mass lost during the RSG phase by comparing the mass at the 
end of the MS phase, 19.5M , with the final mass of the pro- 



genitor, 6M . This gives a total mass loss of 13.5M . The RSG 
life time is 8.75 ■ 10 5 years which implies a mass loss rate of 
1.54 • 10~ 5 M yr _1 . Since we assume that all the mass is lost 
during the RSG phase, we also have to assume that during the 
short duration WR phase the total amount of mass lost to the 
CSM does not increase significantly (see Table[T|i. 

The radii of the forward and reverse shock o f Cas A are 
2:55 + 0:2 and i:58 + 0:i6 (iGotfhelf et al.l l200ll). which at a 
distance of 3.4+j} 3 kpc dReed et all 119951) . correspond to 2.52 
+ 0.2 pc and 1.58 + 0.16 pc respectively. Since the amount of 
mass within the radius of the forward shock has to be ~8M 
(IVinket al.ll 19961) . we can determine the RSG terminal wind ve- 
locity by using the following relation between the mass loss rate, 
wind velocity and forward shock radius: 

8M = / 4nr 2 p(r)dr = / Ixr^dr = , (1) 

in which we assume that the mass loss rate and the velocity 
are constant, and that the density is consistent with a r~ 2 pro- 
file shaped by the free streaming RSG wind. 

When we assume that the density profile, within the cur- 
rent forward shock radius, at the end of the RSG phase is con- 
sistent with an r~ 2 profile shaped by the free streaming RSG 
wind, we have to check if the RSG shell is well outside the cur- 
rent shock rad ius Q For this purpose we used equation (22) from 
IWeaver et a fl dl977l) . which gives the pressure inside the MS hot 
bubble, and the equation for the ram pressure of a stellar wind 
Pram = pv 2 . U sing the parameters f or the MS phase of the stellar 
wind given in Hirs chi et alj (|2004), a typical MS wind velocity 
of 1000 km/s, a ISM density of 10~ 23 g cm~ 3 and the parameters 
for the RSG wind given in Table [1] we obtain a radius for the 
RSG shell of approximately 4 pc. 

Because in our simulations we only consider the inner 3 pc 
of the CSM and the RSG shell is located at approximately 4 pc, 
the assumption of an initial r~ 2 density profile at the end of the 
RSG phase for the inner 3 pc is justified. Since the structure of 
the CSM at the end of the RSG phase can be determined with the 
above mentioned assumptions, it was not necessary to perform a 
hydrodynamical simulation of the CSM during the RSG phase, 
which saved considerable computing time. 

The values for the WR mass loss rate and wind velocity are 
calcu lated with help o f the eq uations (12), (15), (17) and (22) 
from Nugis & Lamers (2000). The parameters needed for the 
calculation of the mass loss rate a nd velocity o f the W R wind 
are taken from the stellar model of IHirschi et alj d2004l) . All the 
values for stellar wind are summarized in Table Q] Since we do 
not know the WR life time, or whether there was a WR phase, 
we treat the WR life time as a free parameter. 

2.2. Numerical method 

All th e simulations pres ented here are done with the ZEUS MP 
code (Ha ves et al.l l2006). This is a three dimensional Newtonian 
magneto-hydrodynamics code, which solves the Euler equations 
on a staggered mesh grid. Magnetism and gravity can also be 
treated by the code, but are not used here. 

In order to separate the CSM from the SN ejecta we use the 
multi-species advection available in ZEUS MP by considering 
two 'species', corresponding to the above mentioned compo- 
nents. The real composition of each of these species is not of 



1 This shell is created by the outflowing RSG wind and a hot MS 
bubble pushing inward. 
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Table 1. Adopted progenitor star wind parameters for our Cas A hydrodynamic simulations: final evolutionary phase of the progen- 
itor star, duration of this phase, corresponding mass loss rate and wind velocity, total amount of mass and kinetic energy lost during 
this phase, and model acronym. 





r[yr] 


M[M />t] 


v[km/s] 


AM[M ] 


A£[10 45 erg] 


Model 


Red Supergiant 


8.75 • 10 5 


1.54 ■ lfr 5 


4.70 


13.5 


5.93 


WR0 


Wolf Rayet 


5 ■ 10 3 


9.72 ■ 1(T 6 


1.716 • 10 3 


0.05 


1.42 ■ 10 3 


WR5 


Wolf Rayet 


10- 10 3 


9.72 ■ 1(T 6 


1.716 • 10 3 


0.10 


2.85 ■ 10 3 


WR10 


Wolf Rayet 


15 • 10 3 


9.72 ■ 1(T 6 


1.716 ■ 10 3 


0.15 


4.27 ■ 10 3 


WR15 



importance since no feedback effects from the species back into 
the hydrodynamic calculations are taken into account. 

Radiative cooling is incl uded using the cooling curve from 
Mac Donald & Bailevl ([198 1), which is valid for a gas of approx- 
imately solar composition, but is applied to all the gas in the 
simulations. Applying this cooling curve to the gas consisting of 
CSM material is reasonable since it has a composition compa- 
rable to the solar composition. Applying it to the gas consisting 
mainly of SN ejecta material is not correct. However, the most 
important cooling in these regions is adiabatic cooling due to 
the expansion of the SN ejecta, and thus the error we make in 
this respect is small. By calculating how much energy we lose 
and assuming that all the gas is optically thin, we can use the 
cooling curve to estimate the radiative energy loss from these 
simulations. 

Our simulations were done in two stages, the first stage was 
that of the interaction of the stellar wind with the CSM, which 
we calculated only once, and the second stage concerned the 
SNR evolution. Simulating the stellar wind was done in the 
same manner as described in lGarcfa-Seg ura et al] d 1996b . In this 
method the innermost radial grid cells are given a density and 
velocity corresponding to the stellar wind at that point in time. 
This is done each time step to simulate the star blowing out its 
stellar wind. We simulated the evolution of the CSM due to the 
stellar wind only once, but used the output of the calculation at 
different evolutionary phases as input for stage 2, the collision of 
the supernova ejecta with the CSM. 

Simulating the supe rnova explosion wa s done in the same 
manner as described in Whalen et al. (2008). This method uses 
the following assumption for the radial density and velocity pro- 
file of the supernova ejecta, which at a given time has a flat inner 
core and a steeply declining outer edge: 

( F ■ r 3 for v < v core 

p(v, = \ F ■ r 3 • (^-T" for v core < v < v max , (2) 
{ Pcsm for v > 

v(r, r) = - for t > . (3) 

Here p is the density, t is the time, v is the velocity and r 
is the radius. F and v core are normalization constants which 
are determined through the assumed SN explosion energy and 
ejecta mass. The maximum velocity (v max ) is set to 3 ■ 10 4 
km/s, which roughly corresponds to the maximum observed 
velocity in core collapse supernovae. The supernova explo- 
sion energy is assumed to be kinetic. The supernova ejecta 
are assumed to be spherically symmetric, so no clumps are 
present in the ejecta. The value for the exponent n is set to 
9 for all simulatio ns, which is usually assumed for core col- 
lapse supernovae dTruelove & McKedll999t iDwarkadai 12005 : 
Chuga i & Chevalier! l2006h ." We calculated four different SNR 
scenarios, corresponding to an explosion inside a CSM shaped 



only by the RSG wind (model WR0), and three explosions in- 
side a CSM shaped by the WR star wind, with WR life times of 
5000 yr (WR5), 10000 yr (WR10), and 15000 yr (WR15). 

For all the results shown, a spherical coordinate system was 
used (i,0,(p), with an assumed symmetry in the ^-direction. The 
number of grid cell in the radial direction was always equal to 
1000. The physical size, in the radial direction, of the grid in the 
pre-supernova evolution simulations was ranging from to 3 pc, 
and the grid cells were equally spaced over this range. 

In the case of the supernova calculations the physical size of 
the radial grid was enlarged each tim e step, which i s also done 
in the same manner as described in Whal en et a D ((2008). For 
each time step the total number of grid cells was redivided over 
the grid, which resulted in an equally spaced grid for the radial 
direction at every time. The number of angular grid points, in 
the ip direction, was equal to 200 for all simulations. The angular 
grid ranged from to | over which the number of grid cells were 
equally spaced. 

Finally, we do not intend to match the observational prop- 
erties of Cas A to the last detail. This would require consider- 
able fine tuning. For example, if we want to match our forward 
shock radius to the observations we could change both the SN 
explosion energy or the SN explosion mass which both affect 
the shock radius. We are more interested in looking at the influ- 
ence of the presence or absence of a WR phase, and its duration, 
on the resulting SNR. 



3. Results 

We have performed only one simulation for the evolution of the 
CSM prior to the supernova explosion. We take the output of 
this simulation at different times and use these different outputs 
as an input for the simulations of the SN ejecta interacting with 
the CSM. 

3.1. Pre-Supernova Evolution 

Figure Q] shows four snapshots of the CSM density structure ap- 
proximately 5000, 10000, 15000, and 19000 years after the start 
of the WR phase. The figure shows that during the WR phase 
not a stable WR shell will form but a fragmented shell in which 
there are several small and larger clumps with a higher density 
than the density in the free streaming RSG wi nd. This clump- 
ing is due to the Vishniac/fhin shell instability ( Vishniac 19831: 
iMac Low & Norman1l993l: lGarcia-Segura & Mac Low 1995b. 

The clumping of the WR shell and its evolution affect the in- 
teraction of the supernova ejecta with the shell. We illustrate this 
in Fig. |2] which shows the cumulative amount of mass for one 
angular grid cell and along all radial cells for the different angles 
of the three models with a WR phase. The relative difference be- 
tween the amount of mass at angles where there is a clump com- 
pared to the angles where there is no clump becomes larger when 
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log(d) [g cm 3 ]: -25.8 -24.7 -23.7 -22,7 -21.6 -20.6 



log d [gem 3 ]: -24.D -23.3 -22.6 -21.9 -21.1 



-24.0 -23.3 -22.6 -21.9 -21.1 




Fig. 1. From left to right the figure shows the CSM density struc- 
ture approximately 5000, 10000, 15000 and 19000 years after 
the start of the WR phase. Looking from the first to the last slice 
(clockwise) the outward moving WR shell can be seen. This 
shell is unstable and will form clumps due to the Vishniac in- 
stability. Outside the WR shell the CSM structure is shaped by 
the wind of the previous RSG phase. 



0.15 




0.05 



Fig. 2. Cumulative amount of mass in one angular grid cell and 
along all radial cells, from to 3 pc, plotted against angle for 
Model WR5, WR10 and WR15. The zero degree angle corre- 
sponds to the left border of the slices of Fig. [1] with the 3 mod- 
els corresponding to the first 3 slices. When the WR phase lasts 
longer with each model, the contrast of the total amount of mass 
at different angles also becomes larger. The peaks in the curves 
correspond to the angles at which clumps are located in the WR 
shell. 



the WR phase lasts longer. This implies that the clump masses 
grow with time. As a result of the clumping the supernova ejecta 
will be able to pass the WR shell easier at certain angles, which 
will affect the shock structure of the resulting SNR. When the 
WR phase lasts longer the clumps in the WR shell will be harder 
to destroy and to accelerate by the SN blast wave. 




X[pc] 



Fig. 3. CSM (blue) and SN ejecta (red) density structure of the 
SNR of Models WR0 (first slice) and WR15 (second slice), 335 
years after the supernova explosion. The blue color corresponds 
to the circumstellar material and red to the supernova ejecta. In 
regions where both material originating from the SN ejecta and 
from the CSM are present, the colors are added and will thus 
become purple. 



3.2. Supernova Remnants 
3.2.1. SNR structure 

Figure [3] shows CSM and SN ejecta density structure of the re- 
sulting SNR for our models WR0 (first slice) and WR15 (second 
slice) 335 years after the supernova explosion. This corresponds 
approximat ely to the age of Cas A assuming an explosion date of 
1671 + 0.9 (iThorstensen et al.ll200lT) . In both slices the forward 
and reverse shock appear as a discontinuous jump in the density. 
The contact discontinuity (CD), which marks the boundary be- 
tween the shocked SN ejecta and shocked CSM, can no longer 
be seen in the second slice, due to the collision of the SN ejecta 
with an irregularly shaped WR shell. In the first slice the fingers 
of the Rayleigh-Taylor instabilities at the CD show the mixing 
of the SN ejecta with the CSM. This instability is caused by the 
high density shocked SN ejecta being decelerated by the low 
density shocked CSM. 

The most important difference between these two models is 
the different structure and corresponding density contrast within 
the SNR. Because the supernova ejecta encountered a smooth, 
spherically symmetric CSM in WR0 the only non-spherical 
component is the one caused by the Rayleigh-Taylor instability 
at the CD. The large density contrasts within the second slice of 
Fig.|3]are caused by the supernova ejecta colliding with a clumpy 
WR shell, instead of a smooth CSM. The high density clumps in 
the WR shell of the progenitor are not completely destroyed by 
the supernova ejecta and can still be seen as high density clumps 
within the SNR. Since they contain a lot of mass compared to 
their surroundings they are hard to accelerate and the supernova 
ejecta moved around them. This is reflected in the structure of 
the shocks, which is somewhat different for angles where there 
is a clump compared to the the angles where there is no clump. 

The distribution of CS mass with respect to the density and 
velocity can be seen in Fig. [4] These figures are mass histograms 
showing the amount of circumstellar mass in a bin correspond- 
ing to a specific range in velocity and density. They were made 
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^wr = o yr 



x WR = 5000 yr 





x WR = 10000 yr 

-21 h 



















log (Mass) [g]: 28.0 28.5 29.1 29.6 30.1 30.6 31.2 31.7 32.2 32.7 



x WR = 15000 yr 




-24.5 - 



-21.5 



-22.5 



-23.5 



-24.5 



2000 4000 

V r [kms" 1 ] 




4000 
V r [km s" 1 ] 



Fig. 4. Circumstellar mass histogram of Model WR0 (top left), WR5 (top right), WR10 (bottom left) and WR15 (bottom right), 
which shows the amount of circumstellar mass in a bin corresponding to a certain range in velocity and density. In all figures the 
high amount of mass at low density and low velocity represents the unshocked CSM. In the upper figures the curved line is an 
artifact from the artificial viscosity, which smears out the forward shock over several radial grid cells. For better comparison we 
have kept the color coding and the axis ranges the same for all subfigures. 



by adding the mass of every grid cell which had a CSM com- 
position in excess of 70%, placing it in the corresponding bin in 
velocity and density. In the upper left mass histogram of Fig. [4] 
there are two regions with a considerable amount of CS mass. 
The first is the unshocked CSM which can be seen as the peak at 
low density and zero velocity. The second is the shocked CSM 
with a velocity exceeding 3500 km/s and a density of approxi- 
mately 10~ 228 g cirT 3 . All the other mass histograms still show 
the unshocked CSM as a high peak in mass at zero velocity and 



low density, but the shocked material is spread out over a much 
larger range in velocities, due to the violent interaction with the 
WR shell. It can also be seen that the range of velocities seen 
in the shocked CSM becomes larger with increasing WR life 
time. This is due to the increasing amount of mass within the 
WR shell for a longer WR life time. Another apparent compo- 
nent in the lower two histograms is the slow moving, highest 
density material which corresponds to the high density clumps 
in the remnant. Model WR5 does not show this feature but the 
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log (Mass) [g]: 28.0 28.4 28.8 29.3 29.7 30.1 30.5 31.0 31.4 31.8 



-21.5 




2000 4000 

V r [kms 1 ] 



6000 



-21.5 



-22 - 



-22.5 



-23.5 



x WR = 15000 yr 




2000 4000 

V r [km s" 1 ] 



Fig. 5. Supernova ejecta mass histograms of Model WR0 (top left), WR5 (top right), WR10 (bottom left) and WR15 (bottom right), 
similar to Figure [4] The freely expanding ejecta can be seen as the horizontal line in all histograms. The rest of the mass in the 
histograms represents the shocked ejecta. 



remainders of the WR shell remain visible even in that model, as 
we will discuss below. 

Figure[5]shows the mass histograms for the SN ejecta for all 
models. The unshocked, freely expanding ejecta are visible as 
a horizontal line in the histograms. Similar to what was seen in 
the shocked CSM distribution, the shocked SN ejecta of Model 
WR0 span a much smaller range in velocities compared to the 
shocked SN ejecta of the models with a WR phase. Also simi- 
lar is that the range of velocities seen in the shocked SN ejecta 
increases with a longer WR life time. 



Figure|6]shows the radiative energy loss of Models WR0 and 
WR15 corresponding to the density plots in Fig. [3] In the fig- 
ure, only the material between the forward and reverse shock is 
visible. The material outside this region does not have a high 
enough temperature to emit enough radiation compared to the 
other components. In the first slice the Rayleigh-Taylor instabil- 
ities dominate, while in the second slice the high density clumps 
are most apparent. WR0 also shows that the RT instabilities are 
composed of both circumstellar material and SN ejecta. In the 
second slice the clumps are by far the most visible feature in the 
SNR, brighter by more than 2 orders of magnitude when com- 
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log'dE,,,,) [erg cm 3 s 21.0 -20.0 -19.0 -18.0 -17.1 -21.0 -20.0 -19.0 -18.0 -17.1 

2.5 - 




Fig. 6. Radiative energy loss of the CSM (blue) and SN ejecta 
(red) components of the SNR Models WR0 (first slice) and 
WR15 (second slice), 335 years after the s upernova explosion. 
The en ergy loss is calculated by using the MacDonald & Bailev 
d!98ll) cooling curve. The second slice clearly shows that due 
to the clumped structure of the WR shell the majority of radi- 
ating material consists of shocked CSM. The fingers due to the 
Rayleigh-Taylor instability in the first slice show emission from 
both shocked CSM and shocked ejecta. 

pared to other material in the SNR. Although we believe that the 
qualitative result of the radiative energy loss is valid, one has 
to be cautious when looking at the results qualitatively. For in- 
stance, non equilibrium ionization and higher metal composition 
would increase line emission and are thus expected to increase 
the cooling rates. 

3.2.2. Shock structure and kinematics 

Figures Q to [10] show the angle-averaged values of the radius 
of the forward shock, the radius of the reverse shock, the ve- 
locity of the forward shock, the velocity of the reverse shock 
and the observed values. The radii of the forward and reverse 
shock were mentioned in Section l2~Tl the v elocity of the forward 
shock is 4990+ 150 km/s dVink et alJl998l:lDeLanev et al.l2004t 
Patn aude & Fe sen 2008) and the velocity of the reve rse shock is 
approximately 2000 ± 400 km/s dMorseet alj|2004l) . The value 
of the reverse shock velocity that is mentioned in the latter paper 
is actually 3000 km/s. However, that is measured in the frame of 
the unshocked ejecta, whereas we are considering the velocity 
of the reverse shock in the observers frame, and should thus use 
2000 km/s, which also comes from that paper. 

In Fig.|7]it can be seen that for a given explosion energy and 
age of Cas A, the longer the WR phase lasts, the larger is the ra- 
dius of the forward shock. As the WR phase lasts longer, the low 
density region blown by the WR wind becomes larger, which en- 
ables the supernova ejecta to expand freely over a longer period 
which increases the forward shock radius. Although the mea- 
surement does not agree with the result of all models, the dif- 
ferences are small compared to the differences in the results for 
the reverse shock which we discuss below. A small change in 
the model parameters, or a difference due to the uncertainties 
within the models themselves, could already change the results 
such that all models could agree with the observations. The rea- 
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Time [yr] 

Fig. 7. Angle averaged radius of the forward shock, for the dif- 
ferent models, plotted against time. The orange dot shows the 
observational value for the forward shock radius at the age of 
Cas A and the line shows the observational uncertainty in ra- 
dius. The uncertainty in the age is too small to be visible in this 
figure. After the initial free expansion phase, the models with a 
WR phase collide with the WR shell at a different time and the 
expansion of the forward shock slows down. 

son for the small variations in shock radius is that all models 
have roughly the same amount of swept up CSM. At the current 
age of Cas A the forward shock velocities of all four models are 
also very similar to one another (Fig. [8]), and lie within 10% of 
each other. 

The comparison of the results for the reverse shock radius 
(Fig. |9]l and the observations tell us more, since there is a clear 
difference between the models. For the models with a WR shell 
there are two important parameters which determine the radius 
of the reverse shock. The first is the radius at which the WR shell 
was located at the time of the supernova explosion, since that 
determines the time frame during which the supernova ejecta 
could expand freely. This is reflected in the different turn off 
times of the initial fast increase in the reverse shock radius for 
the models which contain a WR phase. Model WR0 does not 
expand as fast initially because the supernova ejecta encounter 
a higher density CSM, corresponding to the RSG wind, which 
creates a reverse shock that moves more slowly initially. 

The second important parameter is the amount of mass in the 
WR shell compared to the supernova ejecta mass. Since the latter 
is constant, only the mass in the WR shell is important this case. 
The longer the WR phase lasts, the more mass is accumulated 
in the shell. As a result the collision between this shell and the 
supernova ejecta will be more violent and cause a reverse shock 
which moves inward faster with respect to the forward shock. 
This is represented in the time dependence of the reverse shock 
radii for the different models. The reverse shock radius of WR5 
keeps on increasing after the collision with the WR shell while 
the reverse shock radius of WR15 starts moving inward after 
approximately 300 years. 

This is also visible in Fig. [TO] which shows that the re- 
verse shock velocity of WR15 becomes negative. By comparing 
the models in this figure we can see that for a longer WR life 
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Fig. 8. Velocity of the forward shock, for the different models, 
plotted against time. The observations and corresponding uncer- 
tainty are again shown in orange. Because of the large differ- 
ences in the forward shock velocity during the total evolution of 
the SNR we only plotted a specific range in time. This range 
shows how the values for the forward shock velocity evolve 
around the age at which Cas A is thought to be now. 

time, the velocity of the reverse shock at a given time will de- 
crease, which is what we expect from the arguments given above. 
Together with the reverse shock radius, the velocity shown in 
Fig. [10] shows the largest differences between the four models. 
These differences give us the best tools to constrain the progeni- 
tor of Cas A. 



Fig. 9. Similar to Figure [7] but in this case the reverse shock 
radius is plotted. In orange we show the observational value and 
its uncertainty. Initially all the reverse shock radii are the same 
for the models with a WR shell. However, the WR shell radius 
is different for each model and thus the onset of the decrease in 
expansion, i.e. the collision of the supernova ejecta with the WR 
shell, occurs at different times for WR5, WR10 and WR15, -30, 
~60 and ~ 1 10 years respectively. 



visible in the remnant and are even the dominant emitting com- 
ponent, which we illustrate by showing the radiative energy loss 
of WR5 and WR10 in Fig.[TT] This means that even for a very 



4. Discussion and conclusions 

We have shown the results of 2-D simulations for the hydro- 
dynamic evolution of a SNR in which, for the first time, the 
SN ejecta was separated from the CSM. With the Cas A SNR 
in mind, we investigated the influence of the progenitor wind 
and the interaction of the supernova ejecta with the wind shaped 
CSM on the clumping and shock structure of the SNR. 

From our simulations we can draw the following conclusions 
which can be compared with the observations of Cas A: 

- When a WR phase occurs, the SNR shows the remainders of 
the WR shell. 

- With an increasing WR life time, the reverse shock velocity 
in the observers frame, at the current age of Cas A, decreases. 

- The longer the WR life time, the larger the range of velocities 
of both the shocked SN ejecta and the shocked CSM. 

The clumping that we find in the result of our SNR mod- 
els that include a WR phase, is not comparable to the QSF seen 
in Cas A. In the highest density material we find a somewhat 
larger velocity range, extending to ~ 1000 - ~ 2000 km/s for 
Model WR10 and WR 15, depending on what density range one 
would consider as a clump. This is inconsistent with the obser- 
vations for the QSF. For Model WR5 the clumps are destroyed, 
since no slow moving high density component is visible (Fig. 
[5]). Nevertheless, the remainders of those clumps should still be 



4000 



3000 - 



r~ 2000 
'in 

£ 

> 1000 



- 



-1000 



WR0 
WR5 
WR10 
WR15 

Observations 



280 



300 



320 340 

Time [yr] 



360 



38 



Fig. 10. Similar to [8] but in this case the reverse shock velocity 
is plotted. The models with a longer WR life time have a more 
massive WR shell which causes a higher inward velocity of the 
reverse shock. The model with the highest WR shell mass should 
thus show the largest inward reverse shock velocity in our frame, 
which is consistent with our results. 
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short WR life time of 5000 years remainders of a WR shell will 
be seen in the SNR. 

The clumps in our results are located in between the forward 
and the reverse shock, and their expansion in a 3D picture would 
correspond to a slowly expanding shell. Since all the emitting 
components of our result lie in between the two shocks, they 
should all show a similar distribution when projected onto the 
sky. This is in contrast with what is seen in the observations. Fig. 
3 in lLawrence et al.l (1 19951) shows that the QSF distribution does 
not coincide with the distribution of the other emitting compo- 
nents of the Cas A SNR, and that the distribution of the QSF does 
not have a spher ical shape. Th e same can be seen by comparing 
Fig. 5 and 10 in lFesenl ([2001). This comparison also shows that 
there are some QSF outside the main shell in the southwestern 
part of the remnant, which does not occur in the models. Also, 
the spatial scale of the clumps in our results can range up to 0.2 
pc, which is inconsistent with the measured sizes of the QSFs in 
Cas A. The cooling times of the clumps in our results are of the 
order of hundreds of years (T coo n, lg = % > 100 yr), which does 
not correspond to a typical QSF life time of 25 years found by 
Ivan den Bergh & Kamperl dl985l) . 

To be able to explain the QSF in Cas A we would have to 
invoke another source of dumpiness, since our models are not 
able to explain the QSF. This other source would have to consist 
of clumps with properties similar to the QSF currently seen in 
Cas A. The small size of the QSF in Cas A tells us that they do 
not significantly influence the large scale dynamics of the SNR. 
The effect of the clumps from this other source would thus not 
affect the shock radii and velocities to a large extent and our 
conclusions would remain valid. 

Due to our multi-species approach it can be clearly seen 
that clumps of SN ejecta can be found very close to the for- 
ward shock in models with a WR phase. This is solely due to 
the presence of a clumpy shell. Although it is of no specific sig- 
nificance to our current investigation, there are other remnants in 
which similar features are also observed and explai ned by invok- 
ing cosmic ray acceleration at th e forward shock dWarren et alj 
2005s ICassam-Chenai et afl2 008). While our models and results 
do not apply to these specific SNRs, the fact that a clumpy CS 
shell can have this effect might be of interest to those investiga- 
tions. 

Cosmic ray acce leration is also invoked by 
Patna ude & Fesenl d2008l) . who find that they need the en- 
ergy loss due to this process to be able to explain the shock radii 
and velocities of Cas A. We have not included this process in our 
simulations. Within our models there is a large set of parameters 
to be considered, and we have shown that within this parameter 
space it is possible to find radii and shock velocities which are 
largely consistent with the observed values. This does not mean 
we believe cosmic ray acceleration is not of importance, but 
given the uncertainty in the amount of energy involved in this 
process one can not unequivocally conclude that cosmic rays 
are needed in order to explain the structure and kinematics of 
Cas A. 

The differences in the forward shock radii and velocities of 
our results are not really meaningful, as we could accommodate 
those by small changes in explosion energy and ejected mass. 
However, the reverse shock structure varies considerably in our 
models, and can be used to constrain the progenitor model. The 
radii of the reverse shock in Models WR0 and WR15 lie closest 
to the observed value. From this figure alone one could argue 
that either a WR life time between 10000 and 15000 years or 
no WR phase at all would give us the required shock radius. By 
comparing the reverse shock velocity of WR10 and WR15 with 



log(dE„J [erg cm ' s '), 21.0 -20.0 -19.0 -18.0 -17.1 -21.0 -20.0 -19.0 -18.0 -17.1 




Fig. 11. Radiative energy loss of the SNR of WR5 (first slice) 
and WR10 (second slice), similar to Figure [6] In both slices the 
clumps, which are the remnants of the WR shell, can be seen as 
the strongest emitting component. 



the observations it can be seen that the difference is too large and 
that those models cannot explain the reverse shock velocity seen 
in the observations. The reverse shock velocity of WR0 does not 
match the observed value either, but the difference is smaller. 
If we want to match the reverse shock radius and the reverse 
shock velocity at the same time, a model in which there is no 
WR phase seems most likely. However, recently there has been 
a debate about the reverse shock velocity in the western part of 
Cas A, where the velocity is almost zer o in our frame, and w ould 
be more in line with a long WR phase (Helder & Vink 2008). 

Measurements of th e reverse shock vel ocity are sparse and 
the value mentioned in Mors e et al. I (12004 is rather uncertain. 
Fig. 4.6 of lDeLanevI (2004) shows that th e proper motions of X - 
ray emitting c omponents, which migh t dHelder & Vi nk 2008), 
or might not dPatnaude & Fesenl 120081) be identified with ma- 
terial shocked at the reverse shock, show large differences in 
their expansion rates. In the west of the remnant the expansion 
rate seems to be smaller (-0.1% ~ 0.1% yr -1 ) than that of 
the majority of the material (~ 0.2% yr -1 ). An expansion rate 
of ~ 0.2% yr~' corresponds to a velocity of the X-ray emit- 
ting components of approximately 3000 km/s, given the reverse 
shock radius of 1.58 pc. This would imply a reverse shock ve- 
locity of approxim ately 2500 km/s, wh ich is slightly higher than 
the value given in iMorse et all d2004t) and would be even more 
consistent with Model WR0. 

Our conclusion regarding the visibility of remainders of the 
WR shell in the SNR only holds for the short WR life times 
considered here. If the WR life time was much longer, the in- 
teraction of the WR wind with other parts of the CSM, the 
RSG shell for instance, would alter the structure of the CSM 
and consequently also the structure and visibility of the SNR. 
Nevertheless, this does not change our overall conclusion, since 
we can exclude longer WR life times on the basis of the shock 
structure and kinematics. 

Because the presence of the QSF in Cas A cannot be ex- 
plained with help of the remainders of the WR shell and because 
the remainders of that shell do not correspond to any other ob- 
served component in the SNR, we can exclude the occurrence 
of a WR phase. This strengthens the growing evidence for not 
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having a WR phase during Cas A's progenitor life dSchure et al.l 
2008: lKrause et al.f 2008: Perez-Rend6neiaL||2009|). 
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